Quantum entanglement of particles on a ring with fractional statistics 
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In this paper we investigate the von Neumann entropy in the ground state of one-dimensional 
anyonic systems with the repulsive interaction. Based on the Bethe-ansatz method, the entanglement 
properties for the arbitrary statistical parameter (0 < k < f) are obtained from the one-particle 
reduced density matrix in the full interacting regime. It is shown that the entanglement entropy 
increases with the increase in the interaction strength and statistical parameter. The statistic 
parameter affects the entanglement properties from two aspects: renormalizing of the effective 
interaction strength and introducing an additional anyonic phase. We also evaluate the entanglement 
entropy of hard-core anyons for different statistical parameters in order to clarify solely the effect 
induced by the anyonic phase. 
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I. INTRODUCTION 

In recent years quantum entanglement has attracted 
more and more attentions because it is not only a key 
resource in quantum information theory but also an es- 
sential concept in condensed matter physics P] . It is well 
known that a great deal of properties of quantum many- 
body systems are closely related with the entanglement 
between particles, and the study of which is very im- 
portant for both the system composed of distinguishable 
particles and the system composed of identical particles. 
However, unlike the system of distinguishable particles, 
for which there are various quantities to define and mea- 
sure entanglement, the definition and quantification of 
entanglement between identical particles is still not very 
clear. Fortunately, the von Neumann entropy introduced 
by reduced density matrix is a good quantity to describe 
quantum entanglement for systems consisting of identi- 
cal particles 0,11, i, H, H, 0, Q . The entanglement en- 
tropy of two identical bosons or fermions in abstract wave 
functions has been studied by the Schmidt decomposition 
[1, H, 0, Hi , while the entanglement between two identi- 
cal interacting trapped atoms in a continuous system was 
studied in Refs.[a 0, Si- If shown that the statistical 
properties of identical particles play very important roles 
in their entanglement behaviors. 

As a natural generalization of boson and fcrmion, 
anyon was proposed to describe the particle obeying frac- 
tional statistics and has been a subject of great inter- 
est in the past decades 0, [13, [HI, Hi, [H, [ij, [H. Al- 
though the concept of anyon arises originally from two- 
dimensional systems [H, [H, [3 related to fractional 
quantum Hall effect (FQH) and high-temperature super- 
conductivity d, [13, [HI, the study of ID anyonic model 
has attract ed g reat theoretical interest fisl . [l^ [TtI . [TsI . 
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vated by possible experiments with cold atoms to simu- 
late the creation and manipulation of anyons [U [13, , 
and the possibility of performing topological quantum 
computation (34] , the properties of ID anyons are un- 
der current research focus. Several studies have been 
devoted to ID anyons with a (5-function potential inter- 
action [H, [13, [iaJTij23^^ the limiting cases of 
hard-core anyonr|2l]2i]2il2i, [11 [11, Particularly, 
the ID interacting anyon model is exactly solvable by the 
Bethe-ansatz method as firstly found by Kundu [16]. De- 
spite the intensive studies, the entanglement properties in 
anyonic systems are rarely studied except for the model 
in the hard-core limit . 

In this work, we investigate the entanglement proper- 
ties of a continuous system composed of N anyonic parti- 
cles with repulsive contact interaction on a ring of length 
L by calculating the von Neumann entropy of the single- 
particle reduced density matrix. In general, it is hard 
to calculate the von Neumann entropy of a continuous 
many-body system analytically. So far, most of the stud- 
ies focus on the two- particle system [1, [3, [1| ■ The inte- 
grability of the exactly solvable many-body system pro- 
vides us the possibility to study the entanglement prop- 
erties of a many-body system analytically. Based on the 
Bethe ansatz solution of the interacting anyonic model 
O [13 > evaluate the one-particle reduced density ma- 
trix firstly and then obtain the von Neumann entropy for 
the arbitrary statistical parameter in the full interacting 
regime. 

This paper is organized as follows. In section II, we 
first give a brief introduction to the model and formulate 
the method. In section HI, we first consider the Bose 
limit and focus on the effect of interacting strength on 
the von Neumann entropy. In section IV, we deal with 
the general anyonic case and discuss the effect of sta- 
tistical parameter k on the entanglement properties by 
calculating the von Neumann entropy for different k. A 
brief summary is given in Section V. 
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II. MODELS AND METHODS 

The second quantized Hamiltonian for the one- 
dimensional anyonic system is formulated as 

+ ^ dx^\ (X) (X) (X) (X) , (1) 

in which m is the mass of anyons and gio 
denotes interacting strength between anyons 
[H, [I3|- Here the field operators obey any- 
onic commutation relations 'i>\{xi)'^\{x2) = 
gi.^.(.i-x.)^t^(^2)^t^(^^)^ and «'a(xi)«'^(x2) = 

Sixi ~ X2) + e-*'"'^'(^i""=2)^'^(a;2)^'A(a;i) with 
e{x — y) = 1,-1,0 for a; > y,x < y, and x — y 
respectively. The model ([ij is known to be exactly solv- 
able [10]. The parameter k characterizes the statistical 
property of the anyonic system with k = and k = 1.0 
corresponding to Bose statistics and Fermi statistics 
respectively. The dependence of entanglement properties 
on both the interaction constant gio and statistical 
parameter k (0 < k < 1) will be considered. 

The eigenvalue problem of Hamiltonian ([1]) can be re- 
duced to the quantum mechanical problem of N anyons 
with d interaction 

HiP{xi,...,xn) ^ E^{xi,...,xn) (2) 

with 

^ = -EaT2+2c E Six.-x,), (3) 

i=l « l<i<j<N 

where the natural unit is used and c = mgiu/h? (c > 0) 
is a dimensionless interaction constant. 

In terms of the exact ground-state wavefunction ■(/;, the 
single-particle reduced density matrix is defined as 

/5i =Tr2^3,...,Jv|V')(^l, (4) 

where the trace means to do integrations over all the po- 
sition coordinates except one of them. The single-particle 
entanglement is quantified by the von Neumann entropy 
as 

5 = -Tr(pilog2pi), (5) 

where p is the one-particle reduced density matrix with 
the normalization condition Trp = 1. In coordinate rep- 
resentation, the one-particle reduced density matrix is 
expressed as 

Pi{x,x') 

_ /(f dx2...dxN[lp*{x, X2-, ...XN)ij{x', X2, ■■■,Xn)] 

dxi...dxN\''p{xi,X2, ...,xnW 



Obviously the ([6]) is Hermitian, i.e., we have pl{x,x') = 
Pi{x' , x). The eigen- equation of the one-particle reduced 
density matrix ([6|) is 

/ dx'pi{x,x')(j)^{x') = X^(t)r^{x), (7) 
Jo 

where are the occupation numbers for natural orbitals 
(f>n{x) which form a complete and orthonormal set of 
functions. The one-particle reduced density matrix is di- 
agonal in the basis of natural orbitals and J2^=i ^ ^■ 
In terms of eigenvalues A,, and the eigenfunctions (pnix), 
we can rewrite 

OQ 

pi{x,x') = ^ A^(/)^(a;)0*(x') 
,7=1 

and the von Neumann entropy then reads 

00 

S = -Y,X,\og^\. (8) 

,,=1 

From the above scheme, we can understand the diffi- 
culties which prevent us from studying the entanglement 
properties of a many-body system. First, the calcula- 
tion of the ground-state wavefunction for a many-body 
system is generally difficult except for some exactly solv- 
able systems. Furthermore, even though the exact many- 
body wave function is constructed, the calculation of the 
reduced density matrix for a large system remains a dif- 
ficult task due to the time consuming to calculate multi- 
dimensional integrals. 

III. THE DEPENDENCE OF ENTANGLEMENT 
ON THE INTERACTION STRENGTH 

In this section, we shall focus on the Bose limit and 
study the dependence of entanglement on the interaction 
strength. The effect of fractional statistics will be dis- 
cussed in the next section. In the Bose limit (k = 0), the 
model is reduced to the well-known Lieb-Linger model 
[3^1 . In this case the field operators ^\{x) and ^ a{x) 
obey boson commutation relations and the wavefunc- 
tion il){xi, ...,xm) satisfies exchange symmetry. In the 
scheme of Bethe-ansatz method [35|, the many-particle 
wave function can be formulated as 

■0(a;i, ...,a;jv) = ^e{Xq^ - Xq^_^)...0{Xq^ - Xq^) 

Q 

X(pQ{Xq^,Xq^,...,Xg^), (9) 

where Q labels the region < Xg-^ < Xg^... < Xg^ < 
L, in which qi, q2, qn is one of the permutations of 
1, 2, A'^, Eg sums over all permutations and 9{x — y) 
is the step function. Here ipQ{xq^ ^ Xq2^ ...^Xg^ ) takes the 
Bethe-ansatz type 

V?q(x,i,...,X5„) = X![^PiP2...p« eyiY>{iT.j{kp^Xg^))] (10) 
p 
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FIG. 1: The occupation numbers of the interacting Bose sys- 
tem for different interaction strength c. 



and kp^ is a 



with ^pip2---p]v — ^pY[j<ii^^pi ~ "t 
set of quasi-momentums determined by the Bethe-ansatz 
equations. Here pi,p2, •■•,PAf nieans one of permutations 
of l,2,...,iV, and ep denotes a + (— ) sign associated 
with even (odd) permutation of P. Using the periodical 
boundary condition, we can obtain Bethe-ansatz equa- 
tions [33| whose logarithmic forms are formulated as 



2njn 



N 

E 



2 arctan 



where nj is a set of integers to determine the eigenstates 
and for the ground state nj = (A^+l)/2-j (1 < j < A^). 
The energy of the system is E — X^jLi ^^'^ the total 

momentum is fc = X^^i 

By numerically solving the Bethe-ansatz equations 
(fTTjl . we can obtain the exact ground-state wavefunc- 
tion, and then the one-particle reduced density matrix 
Solving the eigenvalue problem ([7]) numerically, we 
can get a series of Xi and then determine the entangle- 
ment entropy. For simplicity, we shall discuss the many- 
body system with iV = 4 in the following context. In 
Fig. 1, we show the occupation numbers for the inter- 
acting Bose system composed of four identical bosons 
versus different interaction strength c. At c — 0, only 
the lowest nature orbital is occupied which means all the 
bosons condensate to the ground state. The occupation 
number of the lowest natural orbital Ai decreases with 
the increase of the interaction strength, accompanying 
with the increase of occupation numbers of higher natu- 
ral orbitals. Our numerical results for the dependence of 
ground-state entanglement S on the interaction strength 
are shown in Fig. 2. In order to exhibit the change 
with c in a wide range, the logarithm coordinate for c 
is used in this figure. It is shown that the entanglement 
entropy changes monotonically with the change of inter- 
action constant. When there is no interaction between 
bosons (c = 0), no entanglement exists in the system. 




FIG. 2; The entanglement entropy S versus the logarithm 
of the repulsive interaction constant c for a system with 4 
particles. One can see that 5* monotonically and smoothly 
increases from zero to about 1.846 with the growth of c from 
zero to the infinity limit. 



Along with the growth of interaction constant, the entan- 
glement entropy 5* increases slowly in the weakly inter- 
acting regime, and then goes sharply to 1.74. When the 
interaction gets close to the strongly interacting regime, 
the entanglement entropy S slowly approaches to about 
1.846, which is smaller than 2. Our result is accordant 
with [1, 0] J in which we note the entanglement entropy of 
N identical boson-particle system ranges from S = for 
free-boson state to a maximum S in the infinitely repul- 
sive limit which is smaller than S = log2 N. This can be 
understood as follows: when there is no interaction be- 
tween particles, no correlation exists between bosons and 
the occupation number of lowest-energy state is N; while 
in the strong interaction limit (c — > oo), particles will 
be prevented from occupying the same state and higher- 
energy states should be occupied. 



IV. THE DEPENDENCE OF ENTANGLEMENT 
ON THE ANYONIC PARAMETER 

Now we turn to the dependence of ground state en- 
tanglement entropy on the statistics. For anyonic system 
the many-body wave function shall satisfy the general- 
ized symmetry 



1p{- ■■ , Xi, .. Xj , . ..) — 

where the anyonic phase 



''^p{...,Xj,...,Xi,...), (12) 



e{xi-Xk)- ^ e{xj ~ Xk) 



.k=i+l 



for i < j. Considering the symmetry under coordinates 
reflection, we confine k to [0, 1] in the present paper. The 
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FIG. 3: The occupation numbers for different anyonic param- 
eter K with c — 1. 




wavefunction of anyons takes a similar form with that of 
bosons (k = 0) [H, [13 

Q 

X(f)AipQ{Xq^,Xq2,...,Xqj^), (13) 

where (j)A is an additional anyonic phase part 

= exp(-i^ e{Xq^ - XqJ) (14) 

and ipQ{xq-^,Xqr;, ...,Xq,^) has the same form as that of 
Lieb-Liniger Bose gas PH)) except that now we have 

ApiP2...pN = ^pIlj<ii'^Ki ~ *^P3 + c') with 

c' = c/cos(k/2). (15) 

Similarly, the quasi-momenta ki is determined by the 
Bethe ansatz equations 

N k - k 

kjL = 2njTT— ^ 2arctan^^^— ; — (16) 

under the twisted boundary condition ?/'a(0, X2, ■■■,xn) — 
e^"'^^^~^^^AiL,X2: ---jXn)- The Bethe-ansatz equations 
have the same form as that for Lieb-Liniger Boson gas 
(llip if we replace c with the renormalized interaction 
constant c'. According to ([TS]) the effective interaction 
between anyons depends on the statistical parameter k, 
which increases with the increase of k and approaches oo 
in the fermionic limit (k — » 1.0). 

With the same procedure as presented in the above 
section we obtain the one-particle entanglement entropy 
for various c. As a concrete example, in Figure 3 we 
display the change of occupation numbers with differ- 
ent statistical parameters at a fixed interaction strength 
c = 1. The dependence of entanglement entropy on the 



FIG. 4: (Color Online) The entanglement entropy versus the 
statistics parameter k for the anyonic system with c — 1, 
c — 10 and c = 100, respectively. 

statistical parameters k is displayed in Figure 4 for c = 1, 
10 and 100. It is shown that the entanglement entropy 
increases with the increasing of anyonic parameter k and 
interaction strength c. For the case of c = 10, as k = 
the system reduces to Lieb-Liniger gas and the entan- 
glement entropy S is about 0.535; while as k = 1/2, S 
increases to about 1.730 and reaches 2 in the Fermi limit 
with K — 1.0. The similar behaviors are displayed for 
c = 1 and c = 100. For different c the entanglement 
entropy S converges to logj N when k — tt {Smax = 2.0 
here). 

In order to understand why and how k affects the en- 
tanglement entropy, we compare the wavefunction of any- 
onic gas with that of Bose gas. According to (fTS)) and 
([TT|) . it is easy to find that the set of quasi-momenta kj 
depend on c' and the wave function ipQ changes along 
with the change of k. In addition the phase factor (/ja 
includes the statistical parameter obviously. So the influ- 
ence of statistical parameter k on entanglement contains 
two parts: one coming from the phase (f>A and the other 
coming from wavefunction (pq through the renormaliza- 
tion of effective interaction strength c'. 

Finally in order to clarify solely the dependence of en- 
tanglement entropy on the part of anyonic phase 4'a, we 
investigate the hard-core anyons, which can be studied 
using the anyon-fermion mapping method [23l . [23 |. In 
this situation the effective interaction c' = oo and the 
system has the same set of quasi-momenta kj whatever 
the statistical parameter k is. The fractional statistics 
only contributes a phase (jjA in the wavefunction. Ac- 
cording to the numerical result shown in Figure 5, the 
entanglement entropy of hard-core anyonic gas increases 
monotonically when the statistical parameter k changes 
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FIG. 5: The entanglement entropy versus the statistics pa- 
rameter K for the hard-core anyonic gas with = 4. 



entropy in the full interacting regime (0 < c < oo) for 
arbitrary statistical parameter (0 < k < 1.0). Using the 
Bethe-ansatz method, we obtain the exact ground-state 
wavefunction, and thus the one-particle reduced density 
matrix and the von Neumann entropy. In the Bose limit 
(k = 0) the entanglement entropy increases monotoni- 
cally with the increase of the interaction strength and 
approaches a maximum in the hard-core limit. While 
for anyonic system it is shown that the entanglement en- 
tropy increases monotonically both with the increase in 
the interaction strength and statistical parameter. The 
anyonic gas gets to the maximum entanglement entropy 
log2 iV in the Fermi limit {k = 1.0). The statistical 
parameter k affects the entanglement properties of the 
anyonic system by renormalizing the effective interaction 
strength and introducing an additional anyonic phase in 
the wavefunction. The influence of anyonic phase on the 
entanglement is also clarified by evaluating the entangle- 
ment entropy of hard-core anyons for different statistical 
parameters. 



from the Bose limit k = to the Fermi limit k = 1.0. 



V. SUMMARY 

In summary, we have investigated the ground-state en- 
tanglement of the ID anyonic system with repulsive in- 
teraction by calculating the one-particle von Neumann 
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